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INTRODUCTION 


Geomagnetic  ni cr opulsa tions  at  the  surface  of  the  earth  at 
frequencies  well  below  1  Hz  say  conventionally  be  explained  as  being 
the  near  fiel  ’  effects  of  electric  currants  produced  by  decaying 
hydronagnetic  waves  in  the  lower  ionosphere.  The  spatial  variation 
of  the  electromagnetic  field  along  tne  earth’s  surface  will  depend 
not  only  on  the  frequency  of  the  aicropulsations  and  the  conducti¬ 
vity  structure  below  the  surface  bat  also  on  the  spatial  extent  of 
the  ionospheric  currents. 

Using  Fourier  methods,  the  curved,  near  field  wavefronts  asso¬ 
ciated  with  these  geonag:  etic  nicropulsations  nay  be  synthesized 
from  a  spectrum  of  plane  waves  that  includes  botn  uniform  plane 
waves  (planes  of  constant  phase  coincident  with  planes  of  constant 
amplitude)  and  non-unifora  plane  waves  (noc-coincident  planes  of 
constant  phase  and  constant  amplitude) .  Because  of  the  vary  high 
index  of  refraction  of  the  earth  at  frequencies  below  1  ez,  a  uni¬ 
fora  plane  wave  incident  on  a  plane  earth  propagates  essentially 
vertically  into  the  earth  and  the  vertical  field  components  virtu¬ 
ally  vanish.  For  example,  in  a  semi-infinite  conducting  medium  of 
seawater,  the  ratio  of  the  vertical  to  the  horizontal  magnetic  field 
is  less  th3n  10**  7  for  a  frequency  of  10**3  Hz.  Also,  because  of  the 
very  large  free  space  wavelengths  at  frequencies  below  1  Hz,  the 
spatial  variation  of  the  field  along  the  earth’s  surface  is  practi¬ 
cally  negligible  at  all  angles  of  incidence  and  the  surface  is 
essentially  a  plane  of  constant  phase.  In  the  case  of  a  non-unifora 
plane  wave  spectral  component,  the  surface  of  the  earth  is  a  plane 
of  constant  amplitude  and  the  planes  of  constant  phase  in  free  space 
are  noraal  to  the  earth’s  surface.  This  field  geometry  can  result 
in  a  significant  spatial  variation  of  the  electromagnetic  field 
along  the  surface  and  the  spatial  periodicity  of  this  variation  Is 
in  general  unrelated  to  the  free  space  wavelength.  For  a  horizon¬ 
tally  layered  earth  this  horizontal  spatial  variation  must  be  the 
same  for  eJLl  the  layers  in  order  to  satisfy  the  boundary  conditions 
at  each  interface. 

It  is  the  purpose  of  this  report  to  elucidate  the  e.ffect  of  the 
horizontal  spatial  variation  of  the  electromagnetic  field  on 

(1)  the  attenuation  of  both  the  horizontal  and  vertical 
components  of  geomagnetic  micropulsations  ir-  a  two-layer 
medium  with  plane  interfaces 

(2)  the  ratio  of  the  vertical  to  horizontal  magnetic  field 
at  the  top  surface  of  such  a  medium 

(3)  the  surface  impedance  of  a  two-layer  medium. 
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Numerical  results  are  presented  for  a  seawater  layer  over  a 
conducting  earth  bottom. 

Raff  {I960)  and  Price  (1965)  have  addressed  themselves  to  the 
sane  problem;  however,  in  soae  respects  their  treatments  are  not 
as  general.  Raff  besides  choosing  a  rather  unrealistic  spatial 
variation  parameter  does  not  consider  the  ratio  of  the  vertical  to 
horizontal  magnetic  field  components  and  furthermore  presents 
incorrect  numerical  results  for  the  attenuation  of  the  vertical 
magnetic  field.  Price  considers  only  the  case  of  a  non-conducting 
bottom,  which  besides  being  unrealistic,  greatly  exaggerates  the 
effects  of  the  horizontal  spatial  variation  of  the  field. 


FORM  OF  PLANE  WAVES  IN  LAYERED  MEDIA 

A  detailed  treatment  of  electromagnetic  waves  in  stratified 
media  is  available  in  Wait  (1962).  Some  important  results  in  this 
treatment  are  summarized  in  Xraichraan  (1970). 

For  the  purposes  of  the  present  report,  only  the  transverse 
electric  field  (TE)  case  will  be  considered  since  this  case  is 
applicable  for  sources  consisting  of  horizontal  current  sheets  in 
the  ionosphere.  Thus  fsr  a  general,  single  spectral  component 
plane  wave  incident  on  a  horizontally  stratified  conducting  medium 
in  which  the  top  surface  coincides  with  the  X-Y  plane,  the  Z  axis 
points  into  the  layered  medium,  and  the  current  in  the  ionosphere 
flows  parallel  to  the  Y-axis,  the  electric  field  in  the  rath  layer 
is  given  by 

/_  ~umz  .  ,  u mZ\  -jvx 


=  (a  e~UmZ  +  b 
'  m  ,  m 


where  Uj,2  =  V2  +  Ym2#  Ym2  “  ~“2”aem  +  with  the  real  part  of 

Yn  >  0,  and  v  can  take  any  real  value.  A  time  variation  elat  is 
assumed.  In  (1)  the  amplitudes  of  the  downward  and  upward  traveling 
plane  waves  in  the  rath  layer  are  as  and  bm  respectively,  am  is  the 
conductivity  of  the  layer,  and  ym  and  em  are  its  permeability  and 
permittivity  respectively. 

The  magnetic  field  components  in  the  mth  layer  are 
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Since  the  spatial  variation  along  the  X-axis  is  the  same  for 
all  the  layers,  the  propagation  of  the  horizontal  components  of  the 
electromagnetic  field  in  the  layers  may  be  determined  by  using  the 
well-known  transmission  line  analogy  which  is  based  on  the  one-to- 
cne  correspondence  between  reflection  coefficient  at  an  interface 
and  the  impedance  mismatch  ratio. 


THE  ATTENUATION  OF  PLANE  WAVE  FIELDS  IN  A  TWO-LAYER  MEDIUM 

The  equations  ( 1 )  —  ( 3)  for  the  fields  v/hen  combined  with  the 
transmission  line  analogy  readily  give  the  attenuation  of  the  elec¬ 
tric  and  magnetic  fields  in  a  two-layer  medium.  Of  particular 
interest  is  the  ratio  of  a  field  quantity  at  the  bottom  of  the 
first  layer,  z  =  z^,  to  that  field  quantity  at  the  surface,  z  =  0, 
of  the  first  layer.  In  the  event  of  a  stratification  of  more  than 
2  layers,  the  results  below  will  hold  if  the  second  layer  is 
assigned  an  effective  impedance  that  may  be  calculated  from  the 
formulas  in  Kraichman  (1970)  for  the  impedance  properties  of  a 
layered  medium. 


Horizontal  Magnetic  Field  Attenuation 
It  is  desired  to  find  the  ratio  (z^)  /Hlx  (0)  .  From  (2), 


■WV  ■  —  f-i*'"1*1  *  Vi**"11*1) 

3“)^ 


-U1Z1-.  e“3Vx 


»!.«»  .  -O-  (-a,  ♦  .1P12.-2,^» 

DWPj. 


where  the  reflection  coefficient  at  the  interface  between  layers  1 
and  2  is 

n2  ’  ni 

-  _  “  f  f  \ 


n2  *  n. 


and  the  impedances 
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(It  should  be  noted  that  the  reflection  coefficient  for  the  magnetic 
field  is  the  negative  of  p12  defined  in  (6).)  The  first  term  in  (4) 
represents  a  plane  wave  cf  amplitude  aj  traveling  through  the  top 
layer  from  z  =  0  to  z  =  zi  and  the  second  term  is  the  value  of  this 
plane  wave  at  z  =  zj  after  reflection  at  the  interface  between 
medium  1  and  medium  2.  In  (5)  the  first  term  represents  the  value 
of  the  downward  traveling  plane  wave  at  z  =  0  and  the  second  term, 
which  represents  the  value  of  the  upward  traveling  plane  wave  at 
z  =  0,  is  merely  the  first  term  after  an  excursion  through  the  top 
lay^r,  a  reflection  at  the  bottom,  and  a  return  trip  to  z  =  0.  The 
ratio  sought  then  becomes 


Hlx(2l> 

(1  -  P12)e 

Hlx(0) 

1  -  P12e‘ 

e-*lzl 

“2ulzl 


Vertical  Magnetic  Field  Attenuation 


From  (3) 


Elv{zl} 

(1  +  p12)e 

Hlz<°> 
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e~ulzl 

-2*1*1 


using  the  same  reasoning  as  in  (8) . 


RATIO  OF  FIELDS  AT  THE  SURFACE  OF  A  TWO-LAYER  MEDIUM 

It  is  of  interest  to  derive  an  expression  for  the  ratio  of  the 
vertical  to  horizontal  components  of  the  magnetic  field  at  the  top 
surface,  z  *=  0,  of  a  two-layer  medium.  Using  (l)-(3). 


H1Z(°)  = 
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Since  at  z  =  0 
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it  follows  that 
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(13) 


so  that 


Hiz(01 
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(14) 


It  is  also  of  interest  to  derive  an  expression  for  the  surface 
impedance,  T)s#  of  a  two-layer  medium.  This  impedance  is  defined  as 
the  ratio  Ej„ (0) /Hix (0) .  Using  (1)  and  (2)  and  the  value  of  b^/a^ 
in  (13), 


jMMi  /l  +  P12e~2“lZl  \ 

al  V-1  *  p12e‘2“l2V 


(15) 


VALUES  OF  THE  SPATIAL  PROPAGATION  CONSTANT  V 

A  non-periodic  spatial  variation  of  the  field  along  the  sur¬ 
face  of  the  earth  may  be  synthesized  by  a  superposition  of  sinusoi¬ 
dal  components  in  which  the  spatial  propagation  constant  v  takes 
on  all  values  from  0  to  °°.  The  behavior  of  each  sinusoidal  compo¬ 
nent  is  described  by  the  foregoing  equations  in  which  a  spatial 
wavelength  Xg  may  be  associated  with  each  value  of  v  such  that 
v  =  2ir/A3.  The  value  v  -*■  C  or  As  -*•  00  describes  the  field  due  to  a 
current  sheet  of  infinite  width  and  implies  a  normally  incident 
uniform  plane  wave. 

Ionospheric  currents  of  finite  extent  produce  an  aperiodic 
spatial  variation  of  the  field  along  the  earth’s  surface.  The 
effects  of  this  variation  could  in  principle  be  obtained  by  means 
of  Fourier  transforms;  however,  such  a  procedure  is  not  only  complex 
mathematically,  but  also  depends  on  an  accurate  knowledge  of  the 


* 

i 
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aperiodic  surface  variation  which  is  generally  not  available.  It 
is  often  possible  to  get  a  reasonable  approximation  of  the  effects 
by  applying  the  equations  derived  above  to  a  prominent  sinusoidal 
component  of  the  spatial  variation  if  the  spatial  wavelength  of 
such  a  component  can  be  estimated.  Since  in  most  case:,  it  is  im¬ 
practical  to  determine  the  spatial  variation  of  the  fields,  an 
attempt  will  be  made  to  arrive  at  a  range  of  values  for  V  that  will 
represent  spatial  harmonics  that  are  prominent  in  many  geomagnetic 
fluctuations  and  which  can  give  reasonable  estimates  for  the  field 
attenuations  and  ratios  described  earlier. 

Although  the  plane  surface  model  is  used  in  '•.•riving  the 
previous  equations  for  the  two-layer  medium,  *-'*.*  ?. ;t  that  the 

radius  of  the  earth  is  so  much  larger  than  any  depth  of  the  top 
conducting  layer  considered  and  is  also  larger  than  the  skin  depth 
in  the  bottom  layer  lends  some  validity  to  this  approximation. 

This  is  especially  so  for  horizontal  spatial  variations  which  are 
not  completely  global  in  extent.  The  in-phase  ionospheric  currents 
are  most  likely  limited  to  an  earth  quadrant  and  a  prominent  spatial 
wavelength,  As,  along  the  earth’s  surface  is  then  of  the  order  of 
2  quadrants.  Thus,  a  lower  limit  for  V  may  reasonably  be  taken  as 


-7  -1 

v  =  —  r  3  x  10  m 

A 

s 

An  upper  limit  for  v  way  be  estimated  from  the  spatial  field  varia¬ 
tion  along  the  earth's  surface  due  to  an  extremely  localized  source 
such  as  a  line  current.  This  spatial  variation  on  the  surface  of 
a  semi-infinite  conducting  earth  (0  =  1.25  x  10“^  mho/m)  has  been 
calculated  by  Law  and  Fannin  (1961)  for  a  line  current  at  a  fre¬ 
quency  of  approximately  0.05  Hz  located  at  a  height  of  200  km  in 
the  ionosphere.  A  prominent  spatial  frequency  associated  with  field 
variations  along  the  earth's  surface  may  be  estimated  from  the 
graphical  results  presented  by  these  authors.  Thus  the  half  sinus¬ 
oid  that  most  nearly  fits  the  field  variations  in  the  vicinity  of 
the  line  source  has  a  spatial  wavelength  As  =  1.8  x  106  m  and  a 
corresponding  v  =  2tt/As  =  3.6  x  10~"  m"^-.  From  the  foregoing,  the 
values  of  the  spatial  propagation  constant  v  may  reasonably  be  taken 
to  lie  between  the  limits 

3  x  10“?  <  V  <  3.5  x  10"6  (16) 

As  a  demonstration  of  the  practical  value  of  this  approach,  the 
maximum  value  of  | Hlz (0) /Hlx(0) ] ,  which  is  a  sensitive  function  of 
v,  is  seen  from  the  graphical  results  of  Law  and  Fannin  to  be  equal 
to  0.48  and  occurs  at  a  surface  distance  of  400  km  from  the  line 
current  source.  Using  the  above  value  of  v  =  3.5  x  10~®  m~^  in  (14) 
for  the  semi-infinite  medium  yields  0.51  for  this  ratio. 
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NUMERICAL  RESULTS  FOR  THE  TWO-LAYER  MEDIUM 

Numerical  results  are  presented  here  for  a  seawater  layer  of 
conductivity  Uj_  =  4  mho/m  and  a  depth  z^  =  200  meters  over  a  semi¬ 
infinite  earth  bottom  of  conductivities  O2  =  4  x  10“2  mho/m  and  4  x 
10”4  mho/m.  All  layers  are  assumed  to  have  the  permeability  of  free 
space.  Oscillation  frequencies  of  10"*2,  10”3,  and  10“4  Hz  are 
chosen  for  the  ionospheric  current  source.  The  value  of  V  =  10**®  m-^ 
is  used  for  the  spatial  propagation  constant  and  corresponds  to 
approximating  the  spatial  field  variation  by  a  half  wavelength  equal 
to  3  x  10®  meters  or  about  1900  miles.  Such  an  intermediate  value 
for  V  should  be  representative  of  a  goodly  number  of  geomagnetic 
micropulsations . 

Whether  or  not  the  horizontal  spatial  variation  will  affect  the 
horizontal  and  vertical  field  attenuations  given  by  (3)  and  (9) 

spectively  will  depend  on  what  contribution  v  makes  to  the  values 
of  ui  and  U2  which  are  defined  by 

u.  -  (Y^  +  V2)1/2  (17) 

u2  =  (Y22  +  V2)172  (18) 

2  2 

where  Yi  »  jwyo0l  and  Y2  *  j“Vo02*  Table  1  gives  the  ratios 
|v2/y^2|  and  |v2/Y22|  for  various  combinations  of  the  two-layer 
model  parameters  described  above. 

Table  1.  Ratios  jv2/Yj2j  and  Jv2/Y22|  for  Various 

T  awm"  Uni^al  ua  +*V|  Mai  r\  —  6  —,—  1 


10'2  Hz 

10-3  Hz 

10~4  Hz 

V2 

-6 

-5 

-4 

a!  *  4 
mho/m 

3.2  x  10 

3.2  x  10 

3.2  x  10 

°2  = 

4  x  10*2  3.2  x  10‘ 

mho/m 


3.2  x  10 


3.2  x  10' 


°2 

4  x  10~4  3.2  x  10 


0.32 


mho/m 
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It  may  be  seen  from  Table  1  that  the  contribution  of  V  is  quite 
small  for  all  combinations  of  the  parameters  chosen  except  for  the 
frequencies  10“4  Hz  and  1G-3  Hz  and  a  bottom  conductivity  of  4  x 
10~4  mho/m. 

For  the  case  where  the  frequency  is  10“^  Hz  and  02  —  4  x 
10“4  mho/m,  taking  "V  =  10“6  m--1  rather  than  V  =  0  results  in  a  52% 
increase  in  the  horizontal  magnetic  field  attenuation  ratio  given  by 
(8) .  For  the  other  significant  case  where  the  frequency  is  10-3  Hz 
and  O 2  =  4  x  10“4  mho/m,  the  increase  is  only  0.4%.  All  other'com- 
binations  in  Table  1  yield  essentially  the  same  value  for  the  hori¬ 
zontal  magnetic  field  attenuation  ratio  as  that  for  uniform  plane 
wave  incidence  (v  z  0) .  The  vertical  magnetic  field  (or  horizontal 
electric  field)  attenuation  ratio  given  by  (9)  is  virtually  unaf¬ 
fected  even  for  those  cases  in  Table  1  in  which  V  makes  a  signifi¬ 
cant  contribution  to  U2- 

The  racio  of  the  vertical  to  the  horizontal  magnetic  field  at 
the  surface  of  a  two-layer  medium  is  a  rather  sensitive  function  of 
v  since  the  expression  in  (14)  contains  a  factor  that  depends  on  v 
directly.  These  ratios  are  presented  in  Table  2  for  the  two-layer 
model  parameters  described  earlier. 

Table  2.  Ratio  I Hiz (0) /H^x (0) |  for  Various  Two-Layer 


Model  Parameters  with  z^  =  200  meters  and 
v  =  10~6  m-1 


10~2  Hz 

10”3  Hz 

10-4  Hz 

°2  =  4  X 

10-2 

mho/m 

.0091 

.04.; 

.16 

0  =  4  x 

l 

0 

H 

mho/m 

.015 

.13 

.75 

It  is  evident  from  (14)  and  Table  2  that  for  a  given  bottom  conducti¬ 
vity  the  ratio  | Hlz (0) /Hls (0) J  approaches  unity  as  the  frequency 
decreases. 

A  significant  change  in  the  surface  impedance  due  to  the 
horizontal  spatial  variation  occurs  only  for  the  case  in  Table  1  in 
which  the  frequency  is  10“4  Hz  and  the  bottom  conductivity  is  4  x 
10-4  mho/m.  Osing  (15),  the  effect  of  taking  v  =  10“6  m-3-  rather 
than  v  =  0  results  in  a  16%  decrease  in  the  surface  impedance. 
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SUMMARY  AND  CONCLUSIONS 

Spatial  variations  of  geomagnetic  micropulsaticn  fields  along 
the  surface  of  a  horizontally  stratified  conducting  earth  are  the 
result  of  near  field  effects  produced  by  the  ionospheric  current 
sources  of  the  micropulsations.  These  spatial  variations  are  asso¬ 
ciated  with  the  non-uniform  plane  waves  in  the  Fourier  spectrum 
representation  of  the  near  field. 

Assuming  a  two-layer,  horizontally  stratified,  conducting 
medium  and  using  both  a  single  micropulsation  frequency  and  a  single 
spatial  frequency,  expressions  are  derived  for  the  horizontal  and 
vertical  magnetic  field  attenuation  to  the  bottom  of  the  top  layer, 
the  ratio  of  the  vertical  to  horizontal  magnetic  fi.eld  at  the  sur¬ 
face  of  the  two-layer  medium,  and  the  surface  impedance  of  this 
medium.  Limits  are  obtained  for  the  prominent  spatial  wavelengths 
associated  with  geomagnetic  micropulsation  activity.  Using  a  spa¬ 
tial  wavelength  that  is  most  likely  prominent  in  many  micropulsa¬ 
tions,  numerical  results  are  presented  for  a  seawater  layer  over  an 
earth  bottom. 

Assigning  a  value  of  10~6  m_1  for  the  spatial  propagation 
constant  and  the  values  4  x  I0“2  and  4  x  10~4  mI;o/m  for  the  bottom 
conductivity,  results  are  obtained  for  micropulsation  frequencies 
of  10“2;  10“3,  and  10~4  Hz  in  a  seawater  layer  200  meters  deep. 

These  results  show  that  only  for  the  lowest  value  of  bottom  conduc¬ 
tivity  and  micropulsation  frequency  does  the  spatial  variation  along 
the  sea  surface  contribute  significantly  to  the  attenuation  of  the 
horizontal  magnetic  field  to  the  bottom.  The  effect  is  to  increase 
the  attenuation  by  52  percent.  All  other  combinations  of  micropul¬ 
sation  frequencies  and  bottom  conductivities  yield  essentially  the 
same  value  as  that  for  uniform  plane  wave  incidence  (v  *  0).  The 
vertical  magnetic  field  (or  horizontal  electric  field)  attenuation 
to  the  bottom  is  virtually  unaffected  for  all  combinations  of  the 
parameters  above. 

Only  in  the  case  of  the  ratio  of  the  horizontal  to  vertical 
magnetic  field  at  the  surface  of  the  above  two-layer  medium  does 
the  spatial  variation  along  the  surface  have  a  large  influence.  In 
fact,  in  the  absence  of  conductivity  inhomogeneities  other  than  the 
assumed  horizontal  stratification,  it  is  this  spatial  variation 
that  largely  determines  the  value  of  this  ratio  which  for  very  low 
micropulsation  frequencies  and  bottom  conductivities  can  approach 
unity.  It  should  be  noted  that  the  value  of  this  ratio  at  the  sea 
bottom  is  greater  than  its  value  at  the  surface  since  the  vertical 
magnetic  field  attenuation  to  the  bottom  is  less  than  that  for  the 
horizontal  magnetic  field. 

As  is  the  case  for  the  horizontal  magnetic  field  attenuation, 
a  significant  change  in  the  surface  impedance  due  to  the  horizontal 
spatial  variation  occurs  only  for  the  lowest  value  assumed  for  the 
micropulsation  frequency  and  bottom  conductivity  and  results  in  a 
16%  decrease  in  the  surface  impedance. 
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